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Abstract 



We summarise a selection of results on the inviscid limit of the stochas- 
tic Burgers equation emphasising geometric properties of the caustic, 
Maxwell set and Hamilton- Jacobi level surfaces and relating these results 
to a discussion of stochastic turbulence. We show that for small viscosi- 
ties there exists a vortex filament structure near to the Maxwell set. We 
discuss how this vorticity is directly related to the adhesion model for the 
evolution of the early universe and include new explicit formulas for the 
distribution of mass within the shock. 



1 Introduction 

The Burgers equation was first introduced by J. M. Burgers as a model for pres- 
sureles s gas dynamics . It ha s since provided a tool for studying turbulence in 
fluids (jFrisch and Bed . 120011 ). for obtaining detailed a.sympt oti cs for stochastic 



Schrodinger and heat equations (jTruman and Zhad Il996allbl: lElworthv et al 



2001) and has played a part in Arnol'd's work on caustics (lArnol ^dl.ll989l 



1990 



1992) and Maslov's works in semiclassical quantum mechanics (|Maslov and Fedoriuk 
1981) . It has also been used for studying the formation of the early universe 



in th e Zeldovich approximation and also the adhesion model (jArnoFd et al 



Il982l: IShandarin and ZeFdovichl . [l989l) . A detailed explanation of these appli 
cations as well as a com plete history of the Burgers equation can be found in 
(|Bec and KhaninL l2007h . 

In this article we will summarise a selection of results on the inviscid limit of 
the stochastic Burgers equation and outline some applications of these results 
to turbulence and the adhesion model. We begin in Section 2 with a summary 
of results on deterministic Hamilton- Jacobi theory for the heat and Burgers 
equation. 

In Sections 3 to 5 we present so me geometric and analyti c results first devel- 
oped by Davies, Truma n and Zhao ( Davies et al., 2002, 2005) and later extended 
by Truman and Neate (jNeate and Trumanl . 12005, , 2007a ). These results relate 
the geometry of the caustic, Hamilton- Jacobi level surfaces and Maxwell set to 
that of their algebraic pre-images under the inviscid classical mechanical fiow 
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map $t which will be defined in Section 3. In two dimensions these results 
show that a Hamilton-Jacobi level surface, or Maxwell set can only have a cusp 
where their pre-images intersect the pre-caustic and so can only have cusps on 
the caustic. They also allow us to give conditions for the formation of swallow- 
tails on both caustics and level surfaces which in turn have implications for the 
geometry of the Maxwell set. 

We also introduce a reduced (one dimensio nal) action function which was 
developed by Reynolds, Truman and Williams (jTruman et al. L 12003!) under the 



assumption that only singularities of Ak type occur (jArnol' d. 1992). Using 
this, we can find explicit equations for the caustic, level surfaces and Maxwell 
set and their pre-images. In Section 6 we use this to write down an explicit 
stochastic process whose zeros give 'turbulent times' at which cusps on the 
Hamilton-Jacobi level surfaces appear and disappear infinitely rapidly. 

Finally, in Sections 7 and 8, we summarise results showi ng that the fluid has 
non-ze ro vorticity in some neighbourhood of the Maxwell set ( Neate and TrumanL 



2007a() . We show that this vorticity disappears under the assumptions required 
for the adhesion model for the evolution of the early universe and outline a new 
formula for the mass which adheres to the shock (the Maxwell set). 
Notation: Throughout this paper x,Xo,Xt etc will denote vectors (usually in 
R'*). Cartesian coordinates of these will be indicated using a sub/superscript 
where relevant; thus x = (xi, X2, ■ ■ ■ , Xd), xq — {x^, Xq, . . . , Xq) etc. The only 
exception will be in discussions of explicit examples in two and three dimensions 
when we will use (x,?/) and {xq^uq) etc to denote the vectors. 

2 Elements of Hamilton-Jacobi Theory 

We begin by considering a deterministic classical mechanical system consisting 
of a unit mass moving under the influence of a conservative force, —W. This 
system has Hamiltonian, 

H{q,p) = lp^ + Viq), 

where p,q & M''. Let us assume that the system has a given initial velocity field 
VS'o for some function : R'^ ^ R. 

The evolution of this system will be given by the classical mechanical flow 
map, $s : R"* ^ R'' defined by. 



d^ 
ds2 



= -vy($,). 



with initial condition. 



$0 Id, $0 = VS'o, 



where Id denotes the d-dimensional identity map. Thus, if X{s) is a classical 
mechanical path with X{0) = xq, then, 

X{s) = $,(xo), X{0) = V5o(xo). 
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Usually we also demand that X{t) = x for fixed x and t. If Sq and V are twice 
continuously differentiable with bounded second order derivatives, then there 
exists a caustic time tc > 0, such that for all t € {0,tc) the classical mechanical 
flow map is a diffe omorphism. This is a simple co nsequence of the global inverse 
function theorem ( Abraham and Marsden . (l978l) . Therefore we can define, 

xo{x,t) <i>i\x), 

to be the unique pre-image of the point x reached by the path X{s) at time t. 
If we now define. 



S{x,t) := So{xo{x,t)) 



1 ■ 



X\s)~V{X{s)) ds, 



then it can be easily shown that St{x) 
equation, 

dS 

-^+H{x,\/St)^0, St=o{x)^So{x 



S{x,t) satisfies the Hamilton- Jacobi 



(1) 



We now show how the function St can be used to construct a semi-classical 
solution to a corresponding heat equation (jTrumanl . Il977t iTruman and Zhaol 
1996i[l99i). 



Consider the heat equation for u^{x,t) E K where x E K*^ and t > 0, 

,2 



^^l^^u^ + ^,-^V(x)u^, 



with initial condition. 



u^{x,Qi) = exp 



So{x) 



To(x). 



(2) 



(3) 



Let Bs e M'' be a d-dimensional Wiener process on the space JF, P) with 
E {i3(s)i?(i)} = min(s,i). Define an Ito diffusion E M'' and an Ito process 

F/' E «^ by. 



AX^ = -VSt^s{X^) ds + /ids,. 



X^ = X, 

— T 

Iq — X, 



(4) 
(5) 



where < s < t < tc- The time reversal in St^s allows us to effectively 
consider a diffusion process which will reach the point x at time t. Define 
h{s,uj) := ho{Yg^ (uj) , s) where, 

/io(r/,s) ■.^-^l-'vSt^siYn■ 

Since h satisfies the Novikov condition. 



Ep <! exp ( i 



h (s, w) ds I > < oo, 
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where Ep denotes expectation with respect to the measure P, it follows that, 

Ms -.^ cxp h{u,uj)dBu- ^ /i^(u,w)du^, 

is a martingale with respect to J-g — cr{Bs) and P. 

Using the Girsanov theorem, we can now define a new measure P on {Q,J-'), 

dPH = Mt,MdPH, 

and then, 

Bs '■— I h{u,Lj) du + Bs, 
Jo 

is a Brownian motion with respect to P. Therefore, {YfjBs), where Yj^ is 
defined in dSj), forms a weak solution to equation ([4]). That is, 

dY,^ = -VSt-siYI") ds + ndBs, 

and conseqeuntly, 

Ep {fiXI^)} = Ep {/(n^)} = Ep {MJ{Bs)} . (6) 

It follows from the Feynmann-Kac formula that the heat equation ([2]) has a 
solution given by. 



u'^ix, t)^E^\ ToiYn exp ( ~^i-'SoiYn + p^-' / V{Yn ds ) j> , 
and so by equation 



= Ep |To(Xn cxp (^-^i-^So{X>^) + ^ ds 

V5t_.(X,^) dS, - ^ ^* |V5,_,(X,^)p ds) I (7) 

Now, using Ito's formula, 

+/i / V5t_,(X'')dB„ 
Jo 

and so substituting into equation ([7]) for V5f_s(Xf ) di?^ gives, 
u^'(x,t) = e"^Ep|ro(^nexp(^-i^ A5t_,(X,^)ds 
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But St satisfies the Hamilton- Jacobi equation ([T]), and so, by reversing time in 
the diffusion X^, we liave, 



t) = exp 



(8) 



|to(Xo^) exp (^-i 1^ ASt-siXi') ds 

Using the logarithmic Hopf-Cole transformation (lHoDill950l) . 

v^'{x,t) ^ -^^Vlnu^'{x,t), (9) 

the heat equation ([2]) becomes the Burgers equation for velocity field v^{x,t) S 
where /x^ is now the coefficient of viscosity, 



£E- = 5- + K.vk = !^a„..-v,.. 



(10) 



with initial condition. 



v''(x,0) = VS'o(a;) + 0(Ai'). 



We will be particularly interested in the behaviour of for small values 
of fi. In the remainder of this paper we will focus on the discontinuities that 
develop in u'' as /i ^ 0. 

The convergence factor Tq in the initial condition ^ is related to the square 

root of the Burgers fluid mass density , 



(11) 



For t S (0, tc) it can be seen that mass is conserved, 

total mass = J pt{x)dx = J TQ{xQ)dxo — J po{x) dx. 

The next lemma will be key to our treatment of the solution for the Burgers 
equation. 

Lemma 2.1. Consider the above Hamiltonian dynamical system with Hamil- 
tonian H{q,p) and Hamilton- Jacobi function St satisfying, 



dSt^ 
dt 



Hix,VSt) = 0, StMx)=Soix), 



so that, 
Then, 



X{t) = VSt{X{t)), X{0) = WSo{X{0)). 



exp 



ASsiX{s))ds 





dX{0) 


}- 


dX{t) 



where the right hand side is a Jacobian determinant. 
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In particular it follows from Lemma [2.11 for t G {0,tc), that by considering 
an asymptotic expansion of the diffusion Xjf in the solution to the heat equation 



u'^{x,t) = exp ( ] TQ{xo{x,t)) X 



(1 + 0(m')), (12) 



where xo{x,t) is the unique start point of X'^ with, 

i:° = V5.(XO), X? = x. 
Consequently, the Burgers velocity field is given by, 

v" = v>'{x,t) - WSt{x) + 0{fi^). 

3 The stochastic case 

We now consider the behaviour of a Burgers equation with stochastic forcing. 
That is for w^(a;,i) G W^, 

_ + („f . V) «^ = ^At;^ - Vy(x) - eVfct(a;)Tyt, (13) 

with initial condition v^{x,^) — V5o(a;) + 0(/^^), where Wt denotes white noise. 
Using the logarithmic Hopf-Cole transformation the Burgers equation 
becomes the Stratonovich heat equation, 

— = ^Au'' + pr''V{x)u^ + -^kt{x)u^' o Wt, (14) 

with initial condition u^(x, 0) = exp {^^^j^^ Tq{x). 
Now let, 

:=i fx^{s)ds- fv{X{s))ds-e f\s{X{s)) dW^, 
^ Jo Jo Jo 

and select a path X with X(t) = x which minimises ^[X]. This requires, 

dX{s) + VV{X{s)) ds + eVks{X{s)) dWs = 0. 

We then define the stochastic action, A{X{Q),x,t) := inf X{t) = x} . 

Setting, 

A{X{0),x, t) := 5o(X(0)) + A{X{{)), x, t), 

and then minimising A over ^(0), gives ^(0) = V Sq{X{Q)). Moreover, it follows 
that, 

St{x) := inf {A{X{{)),x,t)} , 

X(0) 
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is the minimal solution of the Hamilton- Jacobi equation, 
1, 



dSt+ [ -\VSt\' + Vix)] dt + ekt{x) dWt = 0, St=oix) ^ Soix). 



Following the work of Frcidlin and Wentzell (jFreidlin and Wentzelll . 119981 ), 

— /i^ \nu^(x, t) — !■ St{x), 

as — > 0. This gives the inviscid limit of the minimal entropy solution of the 
Burgers equation as v^{x,t) = VSt{x) ( Dafermo^ 12005 ) . 
Define the classical flow map <l>s : M'^ ^ M*^ by, 



d$, + VVi^s) ds + eVfcs($.) dWs = 0, 



$0 



id, 



V5o. 



Since X{t) = a; it follows that X{s) — ($j ^{x))^ where the pre-image 
Xo{x,t) — ^^^{x) is not necessarily unique. 

Given some regularity and boundedness, the global inverse function theorem 
gives a random caustic time ic('^) such that for < t < tciio), the pre-image, 
XQ{x,t), if it exists, is unique and $t is a random diffeomorphism. Thus, before 
the caustic time v^{x,t) = $f (^^^^{x)) is the inviscid limit of a s olution of the 
Burgers equation with probability one ( Truman and Zhaol . Il996alll998h . 

The method of characteristics suggests that discontinuities in v°{x,t) are 
associated with the non- uniqueness of the real pre-image xo{x, t). In the sitution 
we consider, when this occurs the classical flow map $t focusses an infinitesimal 
volume of points dxo into a zero volume dX{t). 

Definition 3.1. The caustic at time t is defined to be the set, 



Assume that after the caustic time tc{i^) > 0, x has n real pre- images, 

{x} = {a;o(l)(a;, t),xo{2){x, t), . . . , XQ{n){x, t)} , 

where each XQ{i){x,t) e M''. Then the Feynman-Kac formula an d Laplace's 

method in infinite dimensions give for a non-degenerate critical point (jPavies and Truman 
19831 [l984l . 



\x,t) = ^6l,exp 



(15) 



where S}j{x,t) :— Sq {xQ(i){x,t)) + A{xQ{i){x,t),x,t) , and 9i is an asymp- 
totic series in ij? . An asymptotic series in yit^ can also be found for v^{x,t) 
( Truman and Zhaol . [l998). Note that Stix) = min{S'^(x,t) : i = 1, 2, . . . , n}. 

Definition 3.2. The Hamilton- Jacobi level surface is the set, 



m 



{x : 5*0(2;, t) = c for some i} . 



7 



As /i ^ 0, the dominant term in the expansion (fT5|) comes from the min- 
imising xo(i){x,t) which we denote Xo{x,t). Assuming xo{x,t) is unique, we 
obtain the inviscid hmit of the Burgers fluid velocity as the minimal entropy 
v'^{x,t) = $t (ioix^t)). 

If the minimising pre-image xo{x,t) suddenly changes value between two 
pre-images xo{i){x, t) and xo{j){x, t), a jump discontinuity will occur in v^{x, t). 
There are two distinct ways in which the minimiser can change; either two pre- 
images coalesce and disappear (become complex), or the minimiser switches 
between two pre-images at the same action value. The first of these occurs as x 
crosses the caustic. When this results in the minimiser disappearing the caustic 
is said to be cool. The second occurs as x crosses the Maxwell set and again, 
when the minimiser is involved, the Maxwell set is said to be cool. 

Definition 3.3. The Maxwell set is, 

Mt = {x: 3xo,xo € K'' s.t 

X = <i>t(a;o) = $t(io), xq ^ xq and A{xo, x, t) = A{xq,x, t)} . 

We illustrate this in one dimension by considering the integral, 

I{x, t)= G{xo) exp dxo, (16) 

where G € C|J°(M), x e M'' and i = -\/— T. Consider the graph of the phase 
function, F(^^ t-^{xo) — F{xo,x,t), as x crosses the caustic and Maxwell set (see 
Figure [J) . 

As we cross the caustic, the critical point at (a) becomes an inflexion which 
disappears causing Xo{x,t) to jump from (a) to (&). This only causes a jump in 
v^{x, t) when the point of inflexion is the global minimiser of F. As we cross the 
Maxwell set, the critical points at xq and xq move so that F(^^^t){xo) = F(^^^t-^{xo). 
If this pair of critical points also minimise the phase function, then the inviscid 
limit of the solution to the Burgers equation will jump. 





Before 


On Cool part 


Beyond 


Ct 


/-\(a)/ 






Mt 

























Figure 1: Graphs of the phase function as x crosses Ct and Mt. 
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4 The reduced action function 



In this section we will find the phase function F in equation (|16l). We briefl y 
summarise some results of Davies, Truman and Zhao (jPavies et all I2OO2II2OO5I) . 
As before, let the stochastic action be defined as, 



A{xo,po,t) 



X{sfds 



V{X{s)) ds + eks{X{s)) dWs 



where X{s) = X{s,XQ,po) e M"* and for s € [0,t] with xo,po £ R'^, 
dX{s) = -\IV{X{s))ds-e\Iks{X{s))dWs, X{Q) ^ xa, X{Q) 



Po- 



We assume X{s) is unique and let J^g den ote the s igma algebra generated by 
X{u) up to time s. It follows from Kunita (jKunital , Il984l ): 

Lemma 4.1. Assume Sq,V E and kt G C^'^ , W^Vkt Lipschitz with Hes- 
sians V^y, V^fcj and all second derivatives with respect to space variables of V 
and kt bounded. Then for po, possibly xq dependent, 



^ixo,Po,t)^Xit)-^-X^iO), 



1,2,. 



The methods of ([Kolokol^tsov et al.l . 120041 ) guarantee that for small t the 
map Po I— > X{t,xo,pQ) is onto for all xq. Therefore, we can define A(xo,x,t) := 
A{xo,Pa{xo, x,t),t) where Po — Po{xo, x,t) is the random minimiser (which we 
assume to be unique) of A{xo,po, t) when X{t, xq^pq) = x. 

Thus, the stochastic action corresponding to the initial momentum VS'o(xo) 
is A{xq,x, t) := A{xo, x, t) + S'o(a;o). 



Theorem 4.2. is the stochastic flow map then, 

d 



dx: 



-[A{xo,x,t)]^0, a =1,2, 



Using this we can create a one dimensional reduced action function. This is 



done by finding a series of functions Xq{xq, . 



, X, t) for decreasing a 



d,d — 1, . . . , 2 by systematically locally solving the equations. 



dA 

dx\ 



At each stage this eliminates one more coordinate from xq until only a; J remains. 
This gives local reducibility on the assumption that d^A/idx")'^ fo r a = 
2,3, ... ,d and also some mild regularity conditions (ITruman et al. I. I2OO3I) . 

Definition 4.3. The reduced action function is the univariate function, 
f{x,t){xo) •= -4(2^07 ^oi-^o^ 2;, i), . . . , Xq(xq, Xq[-), . . . ,Xq (•), X, t), X, t). 
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The Hamilton- Jacobi level surface is found by eliminating xq between, 



A{xo,x, t) = C, \/xaA{xo, X, t) = 0. 

Alternatively, if we eliminate x to give an expression in xo, we have the pre-level 
surface ^^^H^. Similarly the caustic Ct (and pre-caustic ^^^Ct) are obtained 
by eliminating xq (or x) between, 

/ d^jl \ 
det g{xo,x,t)] =0, Va:oA{xo,x,t) = 0. 

The Maxwell set Mt (and pre-Maxwell set <i>^^M() are obtained by eliminating 
Xq and xo (or x and xq) between the four equations, 

VxoA{xo,x,t) ^ 0, \7j:gAixo,x,t)^0, A{xo,x,t) = A{xQ,x,t) = c. 

The pre-images are calculated algebraically and in the case of the pre-level 
surfaces are not necessarily the topological inverse images. This can be done 
in the free case or when the relevant functions are polynomials in all variables 
which is an implicit assumption in what follows. 

For polynomial A, the eliminations involved with the Hamilton- Jacobi level 
surfaces and caustics are fairly simple to complete using the reduced action 
function with resultants and discrimin ants which can be calculated via Sylvester 
determinants ( van der WaerdenL[l949h . The Maxwell set is more complicated to 



find as eliminating pre-images leads to a surfac e involving both rea l and c omplex 
pre-images termed the "Maxwell-Klein set" ( Neate and Truman . 2005f) . It is 



easier to find the pre-Maxwell set and then use the flow map to parameterise 
the Maxwell set. Parameterising in this manner allows one to restrict the pre- 
image of the Maxwell set to have only real values. In the polynomial case we 
have the following lemma. 

Lemma 4.4. Let denote the polynomial discriminant taken with respect to 
X. The set of all singularities is, 



which factorises as, 



^'(^^U/(x.t)(Ai)-c)) = 0, 



k X Bt{x)^ X Ctixf = 0, 



where Bt — is the equation of the Maxwell-Klein set, Ct ^ is the equation 
of the caustic and k is some non-zero constant. 
The pre-Maxwell set is given by, 

^Ai f /($t(xo),t)(4) " /(ft(xo),t)(Al) ' ^ ^ 



{xi~\^y 

The reduced action function can also be used to ident i fy the cool (singular) 



parts of the Maxwell set and caustic (jNeate and Trumanl . [2007g) . 
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5 Geometric Results 

The results in this section are taken from ( Davies et al. . 2002 : Neate and Trumanl . 
2OO5I 2007a ). Assume that A{xo,x, t) is in space variables with det (^ g^a^^ij ^ 7^ 



0. 

Lemma 5.1. Let denote the stochastic flow map and ^^^Tt and Tt be some 
surfaces where if xq £ ^^^^t then x — <i>t(a;o) G Tt- Then, $t is a dijjerentiable 
map from ^^^^t to Tt with Frechet derivative, 

Let nii{x()), nc{x()) and nM(a;o) denote the normal at xq to the pre-level 
surface, pre-caustic and pre-Maxwell set respectively. Using Lemma |5. II we can 
show that: 

Theorem 5.2. The normal to the pre-level surface is, to within a scalar mul- 
tiplier, given by. 



"HM--(^(^j (^^j Xit,xo,VSo{xo)). 

Theorem 5.3. Assume that a point x on the Maxwell set corresponds to exactly 
two pre-images on the pre-Maxwell set, xo and xo- Then the normal to the pre- 
Maxwell set at xo is, to within a scalar multiplier, given by, 

( d^A \ f d'^A 



^X{t, Xo, V5o(xo)) - X{t, Xo, V5o(io)) 

We now consider the two dimensional case. 

Definition 5.4. Let x = x^-y) = (a;i,a;2)(7) denote a curve where 7 is some 
intrinsic parameter (e.g. arc length) with 7 € (70 — (5, 70 + S) for 70 € R and 
6 > 0. Then the curve is said to have a generalised cusp when 7 = 70 

^(70) =(^(70), ^(70)) =0. 
d7 V d7 d7 / 

It then follows from Theorems 15.21 and 15.31 that: 



Theorem 5.5. Assume that in two dimensions at xq G ^^^H^ the normal 
'ri'Y{_{xo) ^ Q so that the pre-level surface does not have a generalised cusp at xq. 
Then, the level surface can only have a cusp at $4(xo) if^t{xo) G Cf. Moreover, 
if, 

the level surface will have a generalised cusp at x. 
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Theorem 5.6. Assume that in two dimensions at xq G $7 Mt the normal 
nM{xo) so that the pre-Maxwell set does not have a generalised cusp at xq. 
Then, the Maxwell set can only have a cusp at ^t{xo) if^t{xo) G C**- Moreover, 
if, 

the Maxwell set will have a generalised cusp at x. 

These results lead to a range of conclusions relating the geometry of these 
curves. In particular, they allow us to characterise when swallowtails will form 
(a swallowtail perestroika) . The appearance of a swallowtail is related to the 
existence of points with c omplex pre-images which are discussed in detail in 
(|Neate and Trumanl . [2005h . 

Corollary 5.7. Assume that atxo G ^^^Hftl^^^Ct, nii{xo) ^ and nc{xo) ^ 
0. Then at $t(a;o) there is a cusp on the caustic if and only if ^^^H^ touches 
$^^Ct at Xq. Moreover, it follows that, xq G ^^^Mt and that ^^^H^ touches 
$j~^Mt at Xq. Also, at ^tixo), Mt will have a generalised cusp parallel to the 
cusp on Ct- 

Corollary 5.8. Assume that atxo e ^^^MtH^^^Ct, nuixo) 7^ andnc{xo) ^ 
0. Then, there is a cusp on the Maxwell set where it intersects the caustic at 
x = $t(xo) and the pre-Maxwell set touches a pre-level surface ^^^H^ at xq. 
Moreover, if the cusp on the Maxwell set intersects the caustic at a regular point 
of the caustic, then there will be a cusp on the pre-Maxwell set which also meets 
the same pre-level surface ^^^^t '^^ another point xq. 

Corollary 5.9. Assume that at xq G ^^^H^n^^^Ct, nB_{xo) ^ andnc{xo) 7^ 
0. Then at ^tixo) there is a point of swallowtail perestroika on the level surface 
if and only if there is a generalised cusp on the caustic Ct at ^t{xo)- 

The results in this section have natural extensions to three dimensions where 
the cusps are replaced by curves of cusps. We now give some two dimensional 
examples. 

Example 5.10 (The generic cusp). We consider a two dimensional determin- 
istic free example (V = Q, e — Q). In general for such a system the flow map is 
given by, 

^t{xo) ^ xo + iVS'o(xo), 

with derivative map D^t{xo) ~ + tV^S'o(a;o)). The pre-level surface is then 
given by the eikonal equation, 

^|V^o(a;o)P + ^o(a;o) = c, 

where the key identity is, 

V.„ |^|V5o(a;o)l' + ^o(a:o)| - (/ + tV2S'o(a;o))VS'o(xo). 

The generic cusp initial condition, So{xo,yo) = x§j/Oj gives a simple cusped 
caustic (see Figure\^. 
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^t^Ct and ^t^H^ 



Ct and 



Figure 2: The generic cusp caustic (dashed) with three level surfaces (solid line) 

Example 5.11 (The polynomial swallowtail). Let V{x,y) = 0, kt{x,y) = x 
and 5'o(a;o,2/o) = Xq + x^hq. The noisy potential does not affect either the pre- 
caustic or pre-Maxwell set. Consequently at time t the noise will have shifted 
the deterministic caustic and Maxwell set by —e W{u) du in the x direction. 
This point will be returned to in Section 6. 

From Theorem 1 5. 61 the cusps on the Maxwell set correspond to the intersec- 
tions of the pre-curves (points 3 and 6 on Figure\B^. But from Corollary ] 5. 8[ 
the cusps on the Maxwell set also correspond to the cusps on the pre-Maxwell 
set (points 2 and 5 on Figure\M and also Figure\^. The Maxwell set terminates 
when it reaches the cusps on the caustic. These points satisfy the condition for 
a generalised cusp hut, instead of appearing cusped, the curve stops and maps 
back exactly onto itself. At such points the pre-surfaces all touch (Figure^ . 

These two different forms of cusps correspond to very different geometric 
behaviours of the level surfaces. From the definition of a Maxwell set it is clear 
that any point on Mt is a point of self intersection of some level surface. Where 
the Maxwell set stops or cusps corresponds to the disappearance of a point of 
self-intersection on a level surface. There are two distinct ways in which this 
can happen. Firstly, the level surface will have a point of swallowtail perestroika 
when it meets a cusp on the caustic. At such a point only one point of self- 
intersection will disappear, and so there will be only one path of the Maxwell 
set which will terminate at that point. However, when we approach the caustic 
at a regular point, the level surface must have a cusp but not a swallowtail 
perestoika. This corresponds to the collapse of two points of self intersection 
and so two paths of the Maxwell set must approach the point and produce the 
cusp (see Figure\^. 
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\ '5 I 1] 

Pre- curves Curves 

Figure 3: The polynomial swallowtail caustic (dashed) and Maxwell set (solid 
line). 




Cusp on Maxwell set Cusp on caustic 



Figure 4: The caustic (long dash) and Maxwell set (solid line) with the level 
surfaces (short dash) through special points. 




Approaching the caustic Approaching a cusp on the caustic 



Figure 5: The caustic (long dash) Maxwell set (solid line) and level surface 
(short dash). 
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6 Recurrence of stochastic turbulence 

Following (jXruman et all . l2003t iNeate and Trumanl . l2007dl ). the geometric re- 
sults of Section 5 can be used to characterise a sequence of turbulent times. 

Definition 6.1. Real turbulent times are defined to be times t at which there 
exist real points where the pre-level surface <^^^H^ and pre-caustic ^^^Ct touch. 

Real turbulent times correspond to times at which there is a change in the 
number of cusps or cusped curves on the level surface H^. In d-dimensions, as- 
suming <I>t is globally reducible, let f{x,t){xo) denote the reduced action function 
and Xt{X) denote the caustic parameterised using the pre-caustic and flow map. 

Theorem 6.2. The real turbulent times t are given by the zeros of the zeta 
process Ct^ where, 

A satisfies, 

^/(x,(A),t)(Ai) = for a ^ 1,2,..., d, (17) 

and xt (A) is on the cool part of the caustic. 

The term 'real' is used in i jNeate and Trunian . to distinguish this form 

of turbulence from 'complex' turbulence where swallowtail perestroikas occur on 
the caustic. We shall not discuss the details of complex turbulence in this article. 

We now consider the stochastic Burgers equation with white noise forcing 
in d-orthogonal directions, 

_ + („M . v)„P ^ _ (18) 

where W{t) — {Wi{t), W2{t), . . . , Wd{t)) is a d-dimensional Wiener process. 
Proposition 6.3. The stochastic action corresponding to the Burgers equation 

m IS, 

A{x^,x,t) = '"^ + ^(x - xo) ■ f W{s) ds-ex- W{t) 







W{s)Au 



2 



S'o(a;o). 



Lemma 6.4. If xl{X) denotes the random caustic for the stochastic Burgers 
tion fl8]} and x°(A) denotes the deterministic caustic (the e = case) then, 

xliX)=x°iX)-e f W{u)du. 



Using Proposition l6.3l and Lemma [6.41 we can find the zeta process explicitly. 
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Theorem 6.5. In d- dimensions, the zeta process for the stochastic Burgers 
equation Hl^ is, 



t ,2 i-t 



Q^f^^o^^yt)i\i)-ex',i\)-W{t) + e^W{t)- Wis)ds~- iWis^ds-c, 

where f^^ f^i'^i) "is the deterministic reduced action function, Xf{X) is the deter- 
ministic caustic and A must satisfy the stochastic equation, 

Va (/(°.o(,),,)(Ai) - ex?(A) • Wit)) = 0. (19) 

Equation (fT51) shows that the value of A used in the zeta process may be 
either deterministic or random. In the two dimensional case this gives, 

= (v./[L„(,),)(Ai) - eWit)) ■ ^(A), (20) 

which has a deterministic solution for A corresponding to a cusp on the deter- 
ministic caustic. 

Using the law of the iterated logarithm, it is a simple matter to show formally 
that if there is a time t such that Cr — 0, then there will be infinitely many 
zeros of Q in some neighbourhood of r. This suggests that the set of zeros of 
Cf a perfect set which can be rigorously proved in some generality ( Revnold^ , 
2002h . 



The intermittence of turbulence will be demonstrated if we can show that 
there is an unbounded increasing sequence of times at which the zeta process is 
zero. This can be done using an ide a of David W illiams and the Strassen form 
of the law of the iterated logarithm ijTruman et al.. 2003). 

Theorem 6.6. There exists an unbounded increasing sequence of times tn for 
which Yt^ — 0, almost surely, where, 

Yt = Wit)- f Wis)ds-^ f \Wis)\^ds, 
Jo ^ Jo 

and W{t) is a d-dimensional Wiener process. 

Corollary 6.7. Let h{t) = (2ilnln(i))-^. If h{t)^t-^ f^^^o^x) t)(^i) ^ and 

d 

h{i)t^^ ^ x1^{X) 0, then the zeta process is recurrent. 
i=o 



7 A vortex line sheet on the Maxwell set 

We now summarise results of (jNeate and Trumanl . [2007ah which show that for 
compressible flow and small viscosity, with appropriate initial conditions, there 
is a vortex filament structure in the neighbourhood of the cool part of the 
Maxwell set. This result is valid for both deterministic and stochastic cases. 
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We first compute the Burgers fluid velocity on the cool part of the Maxwell 
set where typically for x G Mt the pre-images xo{x,t) and xo{x,t) are well 
behaved functions. Recall from equations (dJ) and that for t < tc{ui) and 
small ^, 



u^{x, t) ^ exp 



Stix) 



P?(x)(1 + 0(a.2)). 



For t > tc{i^), the analagous result for is: 

Lemma 7.1. Let x £ Cool(Mt), so that x = ^t{xo) — ^t{xo) where xq ^ xq 
and A{xq ,x,t) — A{xq ,x,t). Then, 



v^{x,t) - 

where, 

pH^) = TQ{xo{x,t)) 



pUx)VSt{x) + pUx)VStix) 
plix)+pUx) 



0{p' 



dxo , . 



Pti^) = Toixo{x,t)) 



dxo 



and St{x) = A{xo{x, t),x, t), St{x) = ^(2:0(2:, t), x, t). 



We denote by v^{x,t) the leading term for the behaviour of v^{x,t) on the 
Maxwell set and choose orthogonal curvilinear coordinates on Mt denoted by 
(^1,^2)- Let the unit normal in a coordinate patch on Mt be denoted n. 

Theorem 7.2. Ifv'^{x,t) is the leading behaviour of v^{x,t) for x € Mt, then, 



v°{x,t) 



V{St{x)+St{x)) 



dn 



ix) 



dSt 
dn 



pnx)+pn^) 

where ^ denotes the normal derivative {n ■ V) on Mt- 

Definition 7.3. The inviscid limit of the vorticity u'^ is defined to be, 

uj"{x,t) ~ V A v^ix,t), 
where v'^{x,t) denotes the leading behaviour ofv^{x,t) as /i — > 0. 

Since the first term in is on Mt, we obtain: 
Corollary 7.4. For x e Cool(Mt), 



(x) I n 



w° = V A 



— 1 - -^(x) ]n} £ T^Mt. 

.pnx)+p-nx) \9n dn 
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From Corollary 17.41 it follows that oj^ ^Q. Hence, we expect that for small 
viscosity even though inititially there was zero vorticity, once a Maxwell set 
appears, the flow is no longer irrotational. 

The above result shows that for small viscosity a vortex filament structure 
will appear in a neighbourhood of the cool part of the Maxwell set. The limit 
of this vortex filament structure is a sheet of vortex lines on the cool part of 
the Maxwell set. We now give the equation of these limiting vortex lines on the 
Maxwell set in terms of orthogonal coordinates (^1,^2)- 

Theorem 7.5. The limiting vortex lines in a coordinate patch of Mt have equa- 
tions, 

= ciplix) + phx)), 
where c is a real constant and ^ denotes the normal derivative {n ■ V) on Mt- 

This confirms that for our initial conditions, for small viscosity, and for com- 
pressible flow, vortex filaments will inevitably appear in a neighbourhood of the 
cool part of the Maxwell set. Given the rotational effects at work in the universe 
perhaps this suggests that we should consider the Burgers equation with vortic- 
ity from the outset. Kinematical considerations and Galilean invariance suggest 
that the appropriate equation is a Bur gers equation with a vector p otential. We 



hope to discuss this in a future paper (jNeate and Trumanl . 12007 



8 The adhesion model 



The adhesion model for th e formation of the ea rly universe is a refinement of the 
Zeldovich approximation ( ArnolM et al. , 1982f) . In the original adhesion model 
there is no noise and a variational principle is assumed which forces the mass 
to move perpendicular to the cool part of the Maxwell set (often referred to as 
the shock) with the same velocity as the Maxwell set itself. This clearly results 
in mass adheri ng to the cool Maxwell set leading t o an a ccumulation of mass at 
certain points ( Bee and Khaninl . 12007 : Bogaevskv . 2004). A simple calculation 
gives the velocity of the Maxwell set as. 



-{y{St{x)+St{x))] 

and a comparison with Theorem 17.21 and Corollary 17.41 reveals that this adhe- 
sion is sufficient to precisely destroy the vorticity near the Maxwell set. When 
adhesion occurs, an interesting ques tion to consider is how r apidly mass will ac- 
crete on the Maxwell set as ^ (|Bec and Khaninl . |2007|) . This would relate 
directly to the mass involved in the formation of galaxies on the shock. Here 
we prove an inequality for the magnitude of the accumulated mass. 
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For this analysis we need to first mollify the Nelson diffusion process intro- 
duced in Section 2, equation ([4]), 

dX^f = V5„ {X^)du + fidB{u), = x. 

This mollification will remove the discontinuities in the drift caused by the 
Maxwell set and caustic. 

Let < t < T, for some fixed T, and let u e {0,t). Assume that we can 
mollify the minimising Hamilton- Jacobi function Su such that, Su — 5™°'' off 
some thin open set t„ (/x) surrounding the cool Maxwell set M^°°' and the cool 
caustic C^°°' where the Lebesgue measure |tu(^)| = O(^), and that V5™°'' is 
uniformly Lipschitz in space and bounded for u € (0,t). 

We can then define a mollified potential, V^°^^{x), corresponding to this new 
system such that, 

Q cmoU 1 

ou 2 

so that as /i ^ 0, {V^^°^^ — V) is a. surface potential. 

This system then has a corresponding Burgers equation, 

7-) moll 2 

Dt 2 

where this is the physically important velocity field in the limit — > 0. 

Let = /z„ where /x„ is a sequence of real values such that /i„ ^ as 
n ^ oo. Let the diffusion associated with (i = (in after mollification be denoted 
by X" and let, 

A"} = {uj: avoids C^°°' U A/f °°' at all times u € (0, i)} . 

It is important to note that, formally at least, 

P(An = lim 



A- 



exp {^^ Xcc°°i(A-i)(^t"ov('«))dw^ 
X exp ^-A j XMCooi(A-i)(Xt"rov(^i)) 



where. 



Cf°°i(A-i) := {y:d(y,Cf°°i)<A-i}, 
Mf°°HA-^) {y:d(y,Mf°°')<A-i}. 

Here trev denotes time reversal so that, 

We need to assume that for sufficiently large n, P(A") > (5 > 0. 

Now assume that at any point x on the cool Maxwell set, M^°°'^ divides 
space into m{x) parts. Therefore, x has m{x) minimising pre-images XQ{i){x, i). 
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(i = l,2, 

by, 



m) and for each pre-image there is a corresponding S = S{i) given 
VS'„(i)(a;) = X(xo, VS'o(xo), u) 



xo=xo(i){x,t). 

We assume that all VS'(i) are bounded and uniformly Lipschitz in space away 
from the caustic. Note that at a regular point of the Maxwell set m = 2 which is 
the case in which we are interested. In this case Xo{i){x, t) for 1 = 1,2 correspond 
to ma tter arriving on different sides of the cool part of the Maxwell set. Arguing 
as in ([Freidlin and Wentzelll . 119981 ) . the Borel-Cantelli and Gronwall's lemmas 
give: 

Lemma 8.1. For u G (0, t) and n e N (i = 1,2) let, 

dX-^iu) = V SZ{i){Xl\u)) du + ^In dBiu), 

and 

dX?iu)=VSuii)iX?iu))du, 
with X^{t) = X^{t) = X. Then, if J2 fil < oo, for < t < T, 



sup \xnu)-X^{u)\ 

Q<u<t 



as n 



oo 



A7 



1. 



We can now give our results on the distribution of mass. We parameterise 
the cool Maxwell set using the pre-Maxwell set and flow map, 



Cool 



Define xo G i^i 



1 i\jCoq\\ 



and M, 



Cool 



if the classical path from xq to A/, 



Cool 



ids a9° 



for all u < t. It now follows from our results in Section 2 that: 



Theorem 8.2. Let the mass adhering to the Maxwell set in time interval (0, T) 
be m{0,T). Then if, 



mo(0,T) 



dt 
T 



dx^ 



dt 



dxl . 



dx] 



we obtain. 



The volume of M^°°^ 



m(0,r) >mo(0,r). 

is zero and m{0,T) is 0(1). Therefore, the mass 
density on Mp°°^ will be infinite in parts of Mp°°^ in the inviscid limit. As we 
stated earlier, the process of adhesion forces the velocity away from Mt to be 
zero, and consequently the adhering particles carry no vorticity. 

We now conjecture what happens when particles hit C^°°^ first. Taking 
limits, mQ{0, T) is the contribution to the mass in the shock from paths with no 
kinks which would be caused by hitting the caustic before time t. 

It would therefore be reasonable to calculate r7ii(0,r), the contribution of 
mass from paths with one kink caused by a single intersection with the caustic 
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before adhesion occurs. In this case we see that the shock causes a compression 
or decompression of mass. This comes from a generahsed Ito formula for non 
functions. 

Let xq — XQ{t,u, Xq) where Xq G k^{u,t) C M and xq £ K^(u,t) the curviUn- 
ear open set, 



K\t)= y K\u,t). 



If we then integrate over K^{t), and use the generahsed Ito formula for a dis 



continuous function d ue to Elworthy, Truman and Zhao ijElworthv et al 
Feng and Zhaol . l2006l ) . we get, 



2005 



dt 



du 



dxlT^ixo) 



dxo{t,u,xl) 



d{t,u,xl) 



exp {-[V„5„(X(w+)) - V„5„(X(w-))]} , 

where Su{X{u^)) is the action evaluated at the appropriate minimising pre- 
image XQ{i){x,t) just above or below the caustic C„ and V„5 is the normal 
derivative. The final factor in the integrand is the compression/decompression 
term coming from the adhesion process. In one dimension this result can be 
cast as a theorem. In higher dimensions it is only so far a conjecture. 

Clearly, the mass involved in the formation of galaxies in the early universe 
would be given by a sum of terms involving to''(G, T) for r = 0, 1, 2, . . .. As 
discussed at the start of this section, the complete adhesion of all matter to the 
cool part of the Maxwell set precisely destroys the rotation discussed in Section 
7. However, if the adhesion were only partial, this would not totally remove the 
rotation, and could explain the formation of spiral galaxies in the early universe. 
We will discuss such properties and d etailed examples in a forthc oming work on 
the Burgers equation with voriticity (jNeate and Trumanl . [2OO7B). 
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